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Abstract 

The beam quality factor M 2 for the fundamental LP01 mode of a step-index fiber 
is calculated in the presence of gain, in a closed form, as a function of the complex 
generalized V-number. It is shown that the M 2 value of a single-mode gain-guided 
fiber laser can be arbitrary large. The results are important for the interpretation 
of the beam quality measurements in recent experiments on single-mode gain-guided 
fiber lasers. 

1 Introduction 

There has been a growing interest in optical fiber lasers that operate based on the gain- 
guiding index-antiguiding (GG+IAG) principle. Unlike a conventional index-guiding (IG) 
fiber, the core of an IAG fiber has a lower refractive index than the surrounding cladding 
and cannot support IG modes. Siegman [1] has shown that in the presence of sufficient 
gain, an IAG fiber can support confined propagating GG modes, which are normalizable 
in the transverse direction. More interestingly, even for arbitrarily large core diameters, 
GG+IAG fibers can operate in a robust single transverse mode. It is therefore desirable to 
take advantage of the GG+IAG principle and scale up the core size while maintaining the 
single-mode characteristic. The large core size is attractive because it can help mitigate 
the unwanted nonlinear optical effects, raise the optical damage threshold, and increase 
the amplification per unit length of the fiber. 

Several experiments have demonstrated GG+IAG in various fiber laser configura- 
tions [21 [3l HIE]. The reported values of the beam quality factor M 2 in these experiments 
are notably larger than unity, even in fibers that are designed to operate as single-mode. 
This should not be surprising considering the substantial departure of the LP01 profile 



of a GG+IAG fiber from a Gaussian-like beam [BJ, as similarly reported in other uncon- 
ventional optical fibers [7J. However, we show that the values of M 2 for the single-mode 
GG+IAG fibers are substantially larger those measured in the experiments. Therefore, 
the not-so-puzzling larger than unity M 2 measurements in these fibers turn out to be quite 
lower than expected from a pure GG+IAG structure. 

We present a closed form expression for the M 2 parameter of the LP01 mode in the 
presence of gain as a function of the complex generalized fiber V-number. The M 2 param- 
eter has been calculated analytically for the LPmn modes in passive fibers [8]. However, 
we are not aware of a derivation in the presence of gain, which requires integrals of complex 
variables and is mathematically more involved. In this paper, we only present the final 
results and details of the analytical derivation will be reported elsewhere. 



2 Gain-Guided Optical Fibers 

In order to study step-index optical fibers in the presence of gain (GG+IG or GG+IAG), 
it is convenient to use a generalized complex V parameter squared pQ defined as 

V 2 = AN + iG. (1) 

The index and gain parameters AN and G are given by 

AN = (2-Ka/X) 2 2n An, (2) 
G = (2ira/\) 2 (n \/2ir)g, (3) 

where uq is the refractive index of the cladding. riQ + An is the refractive index of the core, 
a is the core radius, g is core power-gain coefficient, and A is the vacuum wavelength. For 
a proper choice of AN and G, the core can support an LP01 guided mode in the form of 

J NJ (ur/a)/J {u), r<a 
E(x,y,z ) = < _ (4) 
I NKq (wr/a)/K (w), r>a 

The parameters u and w are complex and satisfy the following two equation, which can 
be used to determine these parameters, given V . 

uJ\{u)/Jq(u) = wJx{w)/Jq(w), (5) 
u 2 + w 2 = V 2 . (6) 

N is an overall constant to be determined from the normalization condition ([7]) assumed 
throughout this paper, 

JJdS\E(x,y,z)\ 2 = l, (7) 

where J J dS = J J dxdy. In general, the value of V determines the total number of 
confined guided modes in a fiber in the presence of gain, which can be zero or higher. 
The single-mode operating regions of the GG+IG and GG+IAG structures in the AN-G 
Cartesian space are presented in detail by Siegman pQ [6J and will not be repeated here. 
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3 Beam Quality Factor M' 



The beam quality factor M 2 is commonly used in experiments to determine whether a 
fiber laser is single-mode or not. For single-mode fiber lasers, the value of M 2 is near unity, 
while M 2 > 1 is indicative of beam contamination with higher order modes. In order to 
calculate the M 2 , we adopt the variance method, which is mathematically rigorous and 
closely resembles the common experimental procedures for the M 2 measurement [9]. 

Consider an optical beam with the electric field profile E(x, y, z) propagating in the z 
direction. The beam center (x) and the standard deviation of the intensity distribution 



a 2 across the x coordinate are 



(x)(z) 



dS x \E(x, y, z) 



a 2 x (z) = IjdS (x-(x)(z)) 2 \E(x,y,z)\ 2 . 



(8) 
(9) 

Since we only consider cylindrically symmetric optical fibers, the results are identical 
for the y coordinate, and we take the liberty in dropping the x subscript (e.g. M 2 ) 
when convenient. It can be shown that the standard deviation in Eq. \§\ in the paraxial 
approximation, obeys a universal free-space propagation rule of the form 



A 2 

Zo) + B ^ z 



(10) 



zq is the coordinate of the output facet of the fiber which does not necessarily coincide 
with the position of the beam waist. Ref. [8] has shown that 



A 



B 



dS[x — (x)(zq) 



E- 



dx 



c.c. 



dS 



dE 
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+ 
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dS{E 



dE* 
dx 



c.c.) 



(11) 

(12) 



where E = E(x,y,zo) is implied in Eqs. [TT1 and [T2l The position of the beam waist zq^ 
and the beam-quality factor M 2 are given by 



z 0x 

Ml 



Zq - irA/(\B) 



V^Ba 2 (z ) 



A 2 . 



(13) 
(14) 



We note that we slightly differ with Ref. [8] in the sign of the frequency term and also the 
definition of the A term. While the integrals leading the evaluation of the M 2 parameter 
can be evaluated numerically, one is faced with the challenge of reliably truncating the 
infinite integrals. It is difficult to choose a reliable truncation radius, a priori, for a 
desired error. Therefore, a closed-form analytical expression is highly desirable. Using the 
parameters defined in Eqs. [T5l and [TBI 



C^ = (n 2 + ^ 2 r/(n 2 
7 = uJi(u)/J (u), 



★2\n 



U 



(15) 
(16) 
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N, A, B , and <t 2 {zq) can be expressed as 



(4vra 2 iV 



2\-i 



i(Ci° ;u + Ci%).Im( 7 ), 



(17) 



fvra 2 iV 2 l 



A 



<cl,u + cUM^f) - <(h,u - cy - H&,u + Ci%). I7I 2 (is) 



2\"1 



(2vriV 
r2\-l 2 



5 = Re[K 2 C 1 %-^ 2 C 1 %)-7] 



(vra 4 7V 2 )-V 2 (z ) = -Si^u - C 3 1 ^)-Im( 7 ) + (2vra 2 iV 2 ) 1 + 2(C 2 1 ;U + (IJ 
+ 4(C 2 ;U -C 2 %)-|7| 2 



(19) 
(20) 



The contour plots of M as a function of A A" and G are shown in Figs. 1(a) , 1(b) , 1(c) 



The dashed lines in these figures identify the threshold values of the dimensionless gain 
parameter G required to produce a confined and amplifying LP01 mode as described in 



Ref. [lj. The results in Fig. 1(a) related to the GG+IG region, show that M 2 remains very 
close to unity over most of the parameter space. A typical single-mode fiber laser operates 
at 4 < AA" < 6, for which M 2 < 1.04, regardless of the value of the G parameter. This is 
consistent with the common intuition that the beam quality of a single-mode fiber laser 
is very good and does not degrade by pumping, unless higher order modes are excited. 
However, the beam quality factor rapidly increases as the ^-parameter of the fiber gets 



closer to the (dashed) LP01 threshold. This can be seen in Fig. 1(b) which is similar to 



Fig. 1(a) but zoomed in near the origin and the close proximity of the LP01 threshold 
line. 

M 2 



The situation is quite different in the GG+IAG region as shown in Fig. 1(c) 



is quite large over the entire parameter space and becomes exceptionally large near the 
(dashed) LP01 threshold. The large value of the M 2 especially near the threshold is the 
result of the long tail of the beam intensity extending all the way into the cladding region 
over a large portion of the GG+IAG parameter space [6]. The situation closely resembles 
the case of a conventional IG fiber (G = 0), where the tails of the beam extend to the 
cladding in the weakly guiding limit of A — > 0, resulting in large values of M 2 as shown 



in Fig. 1(b) We note that the numerical evaluation of M 2 becomes increasingly difficult 
in the regions of large M 2 and proper care must be taken in the sampling and truncation 
of the beam to get an accurate result. 

In order to explore the impact of the long intensity tail in the cladding on the value 
of M 2 , we introduce a Gaussian apodization function of the form exp(— a _1 r 2 /a 2 ) to 
softly truncate the long intensity tail of beam. We multiply the beam profile of Eq. [J] by 
the apodization function and calculate the M 2 numerically. The results are presented in 
Fig. 1(d) where M 2 is plotted as a function of G for AA" = —100. The dotted line is the 
value of M 2 in the absence of the apodization function. The three solid lines represent 
the calculated values of M 2 for the apodized beam and are marked by different values of 
a = 2,4,8, from the apodization function. Fig. 1(d) clearly shows that the truncation of 



the intensity tail of the GG+IAG beams results in a substantial reduction in the calculated 
value of M 2 . 
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Figure 1: Contour plot of M 2 as a function of AN and G in the a) GG+IG region, b) 
GG+IG region but zoomed in near the origin, and c) GG+IAG region. The dashed lines 
represent the LP01 guiding threshold, d) M 2 as a function of G for AN = —100. The 
solid lines relate to the apodized beam and are marked by different values of a from the 
apodization function. The dotted line is in the absence of the apodization function. 
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The above observations are quite important in relating the M 2 values reported in this 
paper, especially for the GG+IAG region of Fig. 1(c) to the experimental measurements. 
In practice, the extended tail of the GG+IAG beam is truncated at the cladding-jacket 
(or cladding-air) interface. For example, for AN = —100 and G equal to twice the 
LP01 threshold value, the field amplitude drops only by 25% from r = a to r = 2a; for 
AiV = —1000, this drop in amplitude is only 13%. Therefore, the impact of the cladding- 
jacket interface on the tail of beam warrants close attention. 



4 Conclusions 

We have shown that M 2 can be substantially larger than unity even for single-mode fibers 
in the GG+IAG region. The large value of M 2 in the GG+IAG region is the result of 
the long tail of the beam extending all the way into the cladding region, the truncation 
of which can lower M 2 substantially. The reported values of M 2 in several experiments 
on GG+IAG fiber lasers are in the range of 1.05-2.0 [2J H [5] . For example, Ref. [2] 
reports M 2 < 2 for a GG+IAG optical fiber with a 100 lira core diameter and a 250 
fim cladding diameter. Similarly, Ref. [3] reports 1.2 < M 2 < 1.5, where the core and 
cladding diameters are 200 fim and 340 fim, respectively. The cladding diameter is not 
much larger than the core diameter in either experiments. Consequently, the long cladding 
tail of the LP01 beam is truncated at the cladding-jacket interface over a large portion of 
the AN-G parameter space in the GG+IAG region. Therefore, it is not surprising that 
these experiments measure such low values of M 2 . 

We expect that very large values of M 2 will be observed in GG+IAG single-mode 
fibers with sufficiently large cladding to core diameter ratios. Measurement of M 2 can be 
challenging for these fibers since the presence of noise can impact an accurate assessment 
of the tail intensity. The existing experiments [2j [31 HI [5] deviate from the pure GG+IAG 
formalism of Siegman [1] in the sense that the cladding-core diameter ratios are not suf- 
ficiently large. Therefore, an accurate M 2 comparison with these experiments is only 
possible if the presence of the cladding-jacket interface is explicitly taken into account. 
Besides the impact on the value of M 2 , one must examine whether the truncated mode 
remains faithful to its GG+IAG nature or is modified into a mode which is primarily 
index-guided by the cladding-jacket index step. 

The analytical results of this paper have been independently verified by direct nu- 
merical computation. In order to facilitate the implementation of the analytical results 
in this paper for the interested reader, we present the numerical values of the analytical 
solution at V 2 = -5 + 7i. We get: w = 0.79722 + 3.04536i, Na = 0.27742, A = -2.01296, 
Ba 2 = 2.5852, <j 2 x {z Q ) = 0.58852a 2 , and M 2 = 1.42611. 
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